Local moment formation in zinc doped cuprates. 
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We suggest that when zinc is substituted for copper in the copper oxide planes of high Tc 
superconductors, it does not necessarily have a valency of 2+ . Rather, the valency of a zinc impurity 
should be determined by its surrounding medium. In order to study this hypothesis, we examine the 
effect of static impurities inducing diagonal disorder within a one band Hubbard model coupled to a 
localised state. We use this model to discuss the physics of zinc doping in the cuprates. Specifically, 
we discuss the formation of local moments near impurity sites and the modification of the transverse 
spin susceptibility in the antiferromagnetic state. 
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I. INTRODUCTION 



The use of a substitutional impurity such as Zn for Cu. 
within the copper oxide planes of the high Tc niaterialstl 
is considerably helpful for gaining insight into both the 
normal and superconducting properties of these sys- 
tems. Zinc is known to produce a universal rate of Tc 
depression in La.Y,Bi and Tl based high temperature 
superconductorsa, lending support to a universal pairing 
mechanism in these cuprates. It is also known to severely 
affect normal state electronic transport propertiespQ for 
example inducing |Semiconducting behaviour in otherwise 
metallic cuprates.aH The magnetic behaviour of the high 
Tc oxides is also modified in the-|jQresence of zinc. Anti- 
ferromagnetic order is destroyecoQ (albeit not as rapidly 
as superconductivity) and there is a dramatic effect on 
dynamical spin fluctuations. Whilst the underdoped 
cuprates are known to possess a spin pseudo gap, Kakurai 
et.aln have shown that unlike pure YBa2Cu306.6, zinc 
doped YBa2Cu2.9Zno.1O6. 6 has no gap like structure in 
the spin fluctuation spectrum, on the contrary there is 
a significant enhancement in spectral weight at low fre- 
quencies. 

It has been suggested that these normal state effects 
and the observed reduction in T^ are related to the for- 
mation of local jHoments observed in the vicinity of the 
zinc dopants Jjiii These moments can be detected by 
NMR techniquesLj Lj and contribute a .Ci,irje law param- 
agnetism to the magnetic susceptibility.E^llj The purpose 
of this paper is to discuss from a theoretical viewpoint, 
the formation of such local moments and how they relate 
to antiferromagnetism and superconductivity. 

Recently there has been a considerable amount of inter- 
est in local moment formation and its effects on magnetic 
and superconductingiehaviour. Yet there are relatively 
few theoretical worksEj~Ej which attempt to show micro- 
scopically how such local moments may form in a two di- 
mensional strongly correlated fermion system. To date, it 
has been generally believed that zinc substitution results 
in a Zn++ ion situated at a copper site in the supercon- 



ducting planes. This zinc ion is in a closed Sd^*^ config- 
uration and hence the Sd holes of the Cu lattice scatter 
off a static impurity of, from their perspective, very high 
energy. We wish to suggest here that one cannot as- 
sume automatically that the zinc atom is in this valence 
state. The valency of an atom is effected by its surround- 
ing medium. Accordingly, it is not clear that the role 
of the zinc 4s orbital can be immediately discarded from 
the conceptual framework, particularly when 3d fermions 
which are experiencing stroitfe Coulomb repulsion are 
present. In fact, Hao et. aZ.E3 have shown in an elec- 
tronic structure calculation of Lai.85Sro.i5Cui_xZnx04, 
that the partial density of states of the 4s orbital is 
peaked right at the Fermi energy. More ab-initio calcu- 
lations are required to produce a clearer picture on this 
point. 

In this work we consider local moment formation in a 
model which reflects the fact that the zinc Sd electrons 
are tightly bound to the nucleus (absence of holes) whilst 
allowing a coupling between the local 4s orbital and the 
itinerant states of the Cu planar network. To this end 
our model possesses both Hubbard and Anderson type 
qualities. We show that, depending on various coupling 
strengths, local moments are able to form within the 
Hubbard gap and above and below the Hubbard bands, 
with spectral weight on and near to the impurity sites. 
Further, we show that the transverse spin susceptibility 
at the antiferromagnetic wave vector Q = (tt, tt) (the lat- 
tice spacing is set to unity) indeed gains spectral weight 
at low frequencies due to the presence of such local mo- 
ments. 



II. THEORY 

Our model Hamiltonian is the following 

where 
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The creation operator c\^ creates a hole in a 3d copper 
orbital at lattice site i with spin a whilst c g. destroys 
a 3d hole with spin a at lattice site j. The first term 
in H^ represents nearest neighbour hopping on a square 
lattice and the second term is the usual Hubbard inter- 
action between 'id holes at the same site. The last term 
is H^ is simply the 4s orbital energy level where e% is 
the energy difference between the copper 3d and zinc 4s 
hole orbitals. The creation operator s^^ creates a hole 
in a zinc 4s orbital with spin a. The index x labels the 
impurity sites while y labels an impurity site's first near- 
est neighbours. The energy e^ is the difference between 
copper and zinc 3d hole states. Here e^ > whilst ej < 
and ji is the chemical potential of the 3d system. Thus 
H^ represents static impurity scattering of 3d holes off 
the zinc sites whilst H^^ reflects a coupling between the 
localised 4s states and the continuum. 

In two dimensions the Hubbard model itself has proven 
to be rather difficult when it comes to theoretical anal- 
ysis. Consequently, the presence of extra terms in our 
initial Hamiltonian makes our analytical task no easier. 
To this end, we employ thCj^Jiite U slave boson method 
of Kotliar and RuckensteinE3 (KR) which, at the mean 
field level, is in reasonable agreement with. Monte Carlo 
data on the one band 2D Hubbard modeo, and has-iad. 
similar success with the three band Hubbard model.E3"E2l 
We use this technique to generate the antiferromagnetic 
long range order which th*, 2D Hubbard model is believed 
to possess at half fillingB3 We are then able to consider 
the effect of the one body terms H^ and H^^ . While our 
treatment is not self-consistent in the sense the we do not 
use the effect of H^ and H^^ to recalculate the various 
bosonic parameters of the KR method, a Hartree Fock 
self consistent study of a similar model shows that ens* 
rections to this approximation are of higher order only.E2l 

To generate the antiferromagnetic host we introduce 
four auxiliary boson fields representing each possible con- 
figuration of the 3d orbital at a particular site.E3 That is, 
a 3d orbital may be empty, singly occupied with a spin a 
or doubly occupied. The fermionic creation operator is 
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There is a conjugate mapping for the fermionic destruc- 
tion operator. The boson fields e^,p^^ and dj represent 



empty, single with spin a an double occupation of the 
3d orbital at a site i respectively. The denominator in 
Eq. (2) ensures that the mapping becomes trivial in the 
limit t/ — > 0. Unphysical states arising from the en- 
largened Hilbert space are eliminated via the following 
constraints: 



and 
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which imply charge conservation and completeness of 
the bosonic operators respectively. Introducing the con- 
straints into the Hamiltonian H'^ via Lagrange multipli- 
ers we have, 
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where g^^^ = ^ia^*ja- 

To accommodate the antiferromagnetic phase of the 
Hubbard model we split the square lattice into two sub- 
lattices A and B. The symmetry of this phase implies 

In terms of the bosons, 

{p\aPAa) = {PssPBa) = {P^Pa) : 



and 
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We also have A^^ — A^g. = A^.. The partition fmic- 
tion corresponding to 7?° can be calculated exactly in 
the fermion sector and at the saddle point in the boson 
sector (where the bosonic fields become c numbers). The 
self consistent mean field equations can then be solved. 
At half filling the chemical potential n = (A,^ -|- A^)/2. 

At the mean field level, the quasiparticle bandstructure 
is given by 



i?k.a = z^VA2 + g2e(k)2 
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where e(k) is the 2D square lattice tight binding band 
structure, A = (A^ — A J/2 and i/ = ± which refers to 
the upper and lower Hubbard bands respectively. The 
quasiparticle amplitudes on each sub-lattice are: 
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where Okiyo- refers to the A sub-lattice and feki/a to the B 
sub-lattice. 

The Green's function of _ff° may be written as, 
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where Gj^^(k, a;) is the Green's function for the contin- 
uum states and is itself a 2x2 matrix because of the sub- 
lattice structure. The Green's function G°^(w) is due to 
the localised 4s state. Specifically, they are; 
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where Hrj refers to the upper and lower Hubbard bands 
respectively. 

We initially consider the effect of one zinc impurity 
positioned on the A sub-lattice. Formally it is simpler 
to first take into account the effect of H^ . The T- ma- 
trix expression for the continuum Green's function in real 
space is: 
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and of course G^^(w) = Gl'{uj). 

Including now the scattering contributions of H^^ , the 
full continuum Green's function of H in real space is: 

r^ (.)-c^^(.) I ^..'G^.^H^^Gi-(^)^-G^-^,(^) 
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and the full local state Green's function can be written 
as 



g:(^) = 
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where 5](a;) is a self energy given by S(aj) — 

T.yy>t'Gll,{u:)t'. 

The poles oiG1j„{uj) and G%{ijj) give the discrete states 
of the system. They are given by the solution of 
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we recover the 



Note that in the limit t' 
discrete 4s state. 

To determine if the discrete states (which are the solu- 
tions of Eq.([l5|)) are localised in nature we can calculate 
the probability of finding a particle with spin a occu- 
pying one of these states at the impurity site or at its 
nearest neighbours through the residues of the appropri- 
ate Green's functions. For example, for the 4s state we 
use Res[G^(a;); tJ — Edis], for the 3d state at the x site we 
have Res[G^3.£,(w); w = Edis] and for the nearest neigh- 
bour 3d states we use Res[G';^y^{Lj);Lj = Edis]- Here Edis 
is the energy of the discrete state in question. 



III. RESULTS 



A. Local moments 



We have used in our calculations a value oi U = 6At 
which produces a gap of 2A = At and a chemical po- 
tential for the half filled 3d antiferromagnetic system of 
/i = — 3.49t. We have also set the impurity 3d orbital en- 
ergy at e|^ = 9t. This enables us to explore as a function 
of e^ — ^ and t' the nature of the discrete states which 
form in the whole system. 

It is appropriate to first discuss the energy structure of 
the various orbital energy levels and quasiparticle bands. 
In Fig.l we show this relationship. The copper 3d hole 
orbital energy is used as a reference. The two quasiparti- 
cle bands are shown as well as where the Fermi level sits 
and the quantity e^ — fj,. Thus when e^ — ^ = one is 
considering a zinc 4s hole level roughly 3.49t below the 
copper 3d hole level. A large value of |e| — /Li| would be 
favourable for zinc to acquire a valency of 2+ as the 4s 
discrete state would become occupied by holes. 

Solutions to Eq.([l5|) for e^ > and for any values of 
e* — ^ and t' produce two discrete states above the up- 
per Hubbard band (in the hole picture). With increas- 
ing ef they become completely localised to the impurity 
3d orbital. Accordingly they represent localised zinc 3d 
electrons. The appearance of gap states or states below 
the lower band depends on the combination of the above 
parameters. 
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FIG. 1. Energy structure of the hole orbital energies and 
bands (shaded regions). 

Without any coupJJng to the 4s state the system pro- 
duces one gap statalj which has a reasonable amount of 
spectral weight on its nearest neighbours, producing a 
local moment. In the limit e!^ — > cxd the spectral weight 
of the gap state is completely removed from the impurity 
site. In general, for fixed e!^ — 9t, a coupling to the 4s 
state produces one or two extra gap states depending on 
t' and e^ — //. These are the remnants of the 4s localised 
states however some of their spectral weight has been 
shifted to the nearest neighbours by differing amounts 
depending on spin, which in effect produces a local mo- 
ment occupying the 4s state. In Fig. 2 we give an example 
of the graphical form of Eq.(|l5|) for a parameter set of 
t' = 0.25i and ef. - u = -2t. 
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We can see by the intersections with the horizontal 
axis in Fig. 2, that two gap states have formed as have 
two states above the hole upper band. These last two 
states are strongly localised in the zinc 3d orbital and 
show that this orbital is filled with electrons. The up spin 
gap state has amplitude mainly on the zinc neighbours 
however the down spin gap state has amplitude mainly in 
the zinc 4s orbital. Thus we have two local moment con- 
tributions. One situated at the zinc site and one spread 
mainly over the nearest neighbours. Both spins are in 
opposite alignment to the antiferromagnetically aligned 
moments of their respective sub-lattices. 

Both localised gap states are able to act as unpaired 
moments which should contribute a Curie like behaviour 
to the temperature dependent susceptibility. With more 
impurity doping and therefore more localised states, the 
antiferromagnetic correlation length will be reduced lead- 
ing to an eventual loss of antiferromagnetic order. Calcu- 
lations by Sen et alS3 have shown that in the presence of 
local moments, the spin-wave stiffness vanishes beyond a 
certain length scale, destroying long range antiferromag- 
netic order. 
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FIG. 3. Local moment m = n| — n^ as a function of 
(ej — /i)/i for the 4s orbital, the 3d orbital at the zinc site(3d2,) 
and the 3d orbital of one of the nearest neighbours(3dj;). 



FIG. 2. The intersections with the horizontal axis give the 
discrete states of the system. Here we show an example for 
i' = 0.25f , e% ~ fi — — 2t and e^ = 9t.The region between the 
first two dotted vertical lines is the lower band and the region 
between the second pair is the upper band. 



For t' = 0.25t, if e^ — ^ is situated in the gap then an 
up spin gap state and a down spin gap state form with 
large amplitudes in the 4s orbital and thus there is no 
local moment contribution from the 4s states. Accord- 
ingly the zinc atom is nearly inert. An additional up spin 
gap state forms with amplitude mainly on the neighbour- 
ing sites producing a local moment. If however e^ — /i is 



situated within a small energy window between — 3i and 
— 1.5i, then the up spin gap state with amplitude mainly 
on the 4s orbital is lost, leaving only the down spin state 
which in effect leads to a local moment at the zinc site. 
There is still a local moment from the remaining up spin 
gap state whose amplitude is found mainly on neighbour- 
ing sites. If e| — /i is too negative then hole bound states 
are formed below the bands in the 4s orbitals. In Fig. 3 
we show the net moment produced by the discrete states 
as a function of e| — /z in units of t for t' = 0.25t. 

It can be seen that the valence of a zinc impurity can- 
not be simply assumed to be 2+. Our purpose here has 
been to show that the valence of a zinc atom can only 
be considered in the presence of strongly correlated 3d 
fermions. 



B. Spin Fluctuations 

It is instructive to examine the effect of localised states 
on spin fluctuation excitations. Attempts at calculating 
the spin fluctuation spectrum of an itinerant antiferro- 
magnet such as the 2D Hubbard model have generally 
involved a mean field plus RPA type approachEll which 
attempts to incorporate Gaussian fluctuations around 
the antiferromagnetic saddle point. This avenue of at- 
tack has been successful in producing low frequency spin 
wave modes in the large U limit with dispersion which 
agrees with lineajL^pin wave theory of the 2D Heisenberg 
antiferromagnet.t3 Here, we do not attempt to calculate a 
renormalised susceptibility by including fluctuations, but 
instead we calculate the modification of the bare suscep- 
tibility of the antiferromagnetic mean field state due to 
the presence of local moments in an attempt to illustrate 
that the presence of local moments causes an increase in 
spectral weight in the spin fluctuation spectrum, partic- 
ularly at low frequencies. 

Even at this level the full calculation is in general quite 
complicated. Thus, as an example we calculate the con- 
tributions from localised states produced purely by the 
impurity energy e^ . Thus we set t' = and so ignore any 
contributions from localised states formed by coupling of 
the conduction holes to the localised 4s state. 

The transverse spin fluctuations of the rigid antiferro- 
magnetic host can be described by a 2x2 matrix 
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where a and b label the sub-lattices and the elements of 
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where m and n are either a or b. The full spin suscepti- 
bility matrix elements are similarly defined. That is, we 
have 
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For simplicity we assume a small concentration of im- 
purities distributed evenly over the two sub-lattices such 
that multiple scattering from different sites may be safely 
neglected. The diagonal continuum Green's functions 
can be written in momentum space as 
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where n^^ is the concentration of impurities on the m 
sub-lattice. The diagonal components of the susceptibil- 
ity can be calculated by carrying out contour integrations 
in the upper and lower halves of the complex plane, mak- 
ing sure to take into account the poles on the real axis 
due to the discrete states. 

In Fig. 4 we have evaluated at Q = (tt, tt) the imaginary 
part of the localised state contributions to the diagonal 
components of x^(Qi ^) up to first order in impurity con- 
centration oirix^ = Tixf, = 0.05 (full curve). We have cal- 
culated the contributions to the spectral weight of spin 
fluctuations by scattering processes within the bands and 
they are negligible. The dashed curves are the sum of the 
imaginary parts of the diagonal components of x°'^(Q, ^)- 
At this level they represent particle-hole excitations only. 
We have set e!^ = 5t which produces for each impurity one 
gap state and two states above the upper Hubbard band 
in the hole picture. It can be seen that the contribu- 
tion from these localised states is an enhancement in the 
spectral weight at low frequencies. 

The increase in spectral weight at low frequencies 
would be further enhanced if one were to include the off 
diagonal processes which contribute to x'^(Q,w) (which 
are due to inter sub lattice excitations), as well as the 
contribution from any extra localised states produced by 
a coupling of the itinerant states to the zinc 4s states via 
t'. 
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not assume that zinc is in the 2+ state. Rather , its va- 
lency is sensitive to the presence of a strongly correlated 
system. We have shown that local moment formation is 
possible on and near the zinc sites and that they pro- 
duce an enhancement of the low energy spectral weight 
of antiferromagnetic spin fluctuations. 
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FIG. 4. Imaginary part of the transverse spin susceptibility 
at Q = (tTjTt). Here e!^ = 5i. The dashed line is the sum of 
the diagonal components of Im x'^'^(Qi '^) whilst the solid line 
represents the contributions from the discrete states. 

Indeed if this increase in spectral weight is also pro- 
duced away from half filling then this phenomena may go 
some way in explaining how the spin pseudo gap in the 
underdoped xompounds vanishes with zinc doping. As 
Tallon et alr3 have pointed out, for a small amount of 
zinc doping the pseudo gap is suppressed near the vicin- 
ity of the dopants and may be completely suppressed 
when the mean spacing between zinc atoms falls below 
the pseudo gap correlation length. 

In this picture superconductivity would also be extin- 
guished when the pair correlation length falls below the 
mean spacing between zinc atoms. Pair breaking could 
occur via the Abrikosov-Gorkov mechanism due to ef- 
fective magnetic scattering, however there is evidence to 
suggest that this effect is too small to account for the 
rapid reduction of Tc.c3 If the order parameter has dj.2_y2 
symmetry then in addition to the this, strong potential 
scattering from the zinc impurities may be enough to 
suppress Tc at a rate seen by experiment. 

Alternatively, if in the underdoped materials, the for- 
mation of a spin pseudo gap is related to the formation 
of local singlet pairs, then suppression of this gap due to 
an increase in spin excitation spectral weight may also 
be a possible explanation for the reduction of Tc- 



IV. CONCLUSIONS 

We have attempted to show that when substituting 
zinc for copper within the copper oxide planes one can- 
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